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Abstract. Assuming that the interchange convection of 
magnetic flux elements is the physical cause for the ex- 
istence of filamentary penumbrae in sunspots, we inves- 
tigate the behavior of an individual fibril embedded in 
the deep penumbra. The fibril is approximated by a thin 
magnetic flux tube which evolves dynamically in the en- 
vironment given by the global magnetostatic model of a 
sunspot. 

Our simulation shows that the flux tube, initially posi- 
tioned at the penumbra-quiet Sun boundary in the model, 
will rise through its deep penumbra developing a flow 
along the tube that points upward beneath the photo- 
sphere, and radially outward above the photosphere. Our 
results suggest that a bright filament may be formed by 
an extended tail of a penumbral grain. Such filaments 
are optically thick, hotter than the surroundings, and ele- 
vated above a darker background. An upflow in penumbral 
grains bends horizontally outwards above the photosphere 
and gradually cools down due to radiative losses leading 
to a tail that gradually darkens. The plasma flow inside 
the flux tube then becomes transparent and the tube con- 
stitutes a thin elevated flow channel, that can reproduce 
the observed features of the Evershed effect. We present 
also a new acceleration mechanism for the Evershed flow. 
It is demonstrated that a local surplus of gas pressure 
develops inside the tube as it rises through the specific 
(superadiabatic and magnetized) penumbral background. 
The resulting gradient of the gas pressure can drive the 
flow along the tube. 
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1. Introduction 

A sunspot penumbra seen in high resolution images reveals 
a small scale fine structure (for a review, see Muller |1992D . 



At a resolution better than 0.5", penumbral grains and 
radially elongated bright and dark filaments become visi- 
ble. There is some observational evidence that bright fila- 
ments consist o f a few penum bral g rains which are radially 
aligned (Muller 1973b ; Soltau 1982 ). Penumbral grains ex- 
hibit proper motions and migrate inwards, i.e. toward the 
umbra , with apparent velocities of 300-500 m s _1 (Muller 
1973a ). Muller describes penumbral grains as comet-like 



structures, having a bright coma and a somewhat dim- 
mer tail, which is always directed radially outwards, i.e. 
away from sunspot center (see also Fig. 2 in Tritschler 
et al. ( 1997 )). He states that penumbral grains have a 
diameter of 0.35". However, this value is close to the spa- 
tial resolution limit that was achieved in his observations. 



Grossmann-Doerth & Schmidt (1981) find on the basis 
of a statistical analysis of penumbral fine structure, that 
bright and dark filaments have to be smaller than 0.55". 
Using the 1.5m MacMath telescope, Stachnik et al. (1983) 
find spatial scales of 0.11" within the penumbra. These 
results rely on power spectra which were obtained using 
speckle interferometry. Thus, since the spatial dimensions 
of bright and dark filaments are uncertain, it is impos- 
sible to assign reliable values for corrected photospheric 
intensities and effective temperatures of the penumbral 
fine structure. 



A still not fully understood penumbral phenomenon 
was discovered by Evershed ( 1909j ): photospheric absorp- 
tion lines show a shift of the line core that turns out to be 
related to a line asymmetry when dark penumbral struc- 
tures are observed (Beckers & Schroter 



1969; Wiehr 



Degenhardt 1994). In recent years evidence has accumu- 
lated that the Evershed effect must be due to a magnetized 
plasma flow (Solanki et al. 1994) occurring in thin elevated 
horizontal channels containing outward flows (Rimmclc 



1995). The observed asymmetry of the intensity profile can 
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then be explained as the superposition of an unresolved 
Dopplcr shifted component and an unshifted main compo- 
nent. Here it is assumed that an individual flow channel 
cannot be spatially resolved as suggested by many ob- 
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servers (Bumba 1960; Holmes 1963; Schroter 1965; Stell- 
macher & Wiehr |l980| ; Wiehr |1995| ; Balthasar et al. |l997| ). 

The siphon flow mechanism has been proposed to ex- 
plain the Evershed effect by Meyer & Schmidt ( |l968[ ). It 
was studied numerically for stationary flux tubes, which 
arch from the inner penumbra to the outer penumbral 
boundar y, by numerous authors (e.g. , Tho mas 1988 ; De- 
genhardt 1991; Thomas & Montesinos 1993). Hitherto, the 
stationary siphon flow model is the most broadly accepted 
explanation for the Evershed effect. For these models it is 
essential to assume strong magnetic flux concentrations 
of opposite polarity near the outer edge of the penumbra 
that can create a gas pressure gradient along the arched 
flux tube in order to accelerate the plasma. Here, we want 
to present for the first time a dynamical model in which 
the plasma is accelerated locally in the penumbra. In our 
model, the gas pressure gradient that accelerates the flow 
is naturally caused by the evolution of the flux tube. 

The subphotospheric structure of sunspot penumbrae 
cannot be observed directly. It is deduced by means of 
numerical models for the overall structure of a spot (see 



the review of Jahn 1992); models with characteristics at 
the surface which are concurrent with observations yield 
some information about the spot's properties in deeper 



layers (e.g., Jahn [1989| ; Jahn & Schmidt |1994[ ). Both mod- 
els and observations of the magnetic field strength and its 
inclination within the penumbra (e.g., Beckers & Schroter 

Title et al. |1993| ; Lites et al. 



196E; Schmidt et al. 1992 



1993) indicate that the penumbra cannot be a surface 



phenomenon. Instead, it has to be deep (Schmidt 1987 



1991) 



The averaged photospheric heat flux of the penum- 
bra, F pu « 0.75Fq, is significantly higher than the corre- 
sponding umbral heat flux, F u ss 0.23Fq. Jahn & Schmidt 
( |1994 hereafter JS) proposed a concept of interchange 
convection of magnetic flux tubes to account for the sur- 
plus brightness of the penumbra. Inspired by observations 
of bright and dark filaments they conjectured that the 
magnetic field within the penumbra consists of an en- 
semble of magnetic flux tubes. They surmised that in- 
clined magnetic flux tubes evolve dynamically within the 
penumbra and are able to transport heat more efficiently 
than magnctoconvection in the umbra which contains an 
almost vertical magnetic field. The onset of interchange 
convection should depend on the slope of the magnetic 
field lines: Imagine a perturbation acting on a bundle of 
magnetic field lines being in magnetohydrostatic equilib- 
rium. Since the component of the buoyancy forces acting 
perpendicular to the tube increases as the flux tube be- 
comes more horizontal, there may exist an angle for which 
buoyancy forces perpendicular to the tube resulting from 
a perturbation overcome restoring magnetic forces. Then 
interchange convecti on of magnetic flux tubes sets in (see 
also Rucklidge et al. 1995 ). 

JS constructed the tripartite model for a sunspot 
which describes self-consistently the three stratifications 




radius in Mm 

Fig. 1. Sunspot model in a meridional cut. Above 
Zbp, the magnetopause transmits heat which is distributed 
horizontally by the interchange of magnetic flux tubes. 
Both the peripatopause and the lower part of the magne- 
topause do not transmit energy (see also Fig. 2 in JS). 



corresponding to the umbra, the penumbra, and the quiet 
sun. These three stratifications are in mechanical equilib- 
rium, such that different gas pressures between the ad- 
jacent stratifications are balanced by the magnetic pres- 
sure. The corresponding magnetic fields are generated by 
two current sheets, namely the magnetopause, which sep- 
arates the quiet sun and the penumbra, and the peri- 
patopause, which separates the penumbra and the umbra. 
Fig. [l] shows a schematic drawing of these current sheets. 

All the stratifications are in hydrostatic equilibrium 
and the heat transport (with varying efficiency) is de- 
scribed by means of the mixing length theory. It is as- 
sumed that the peripatopause thermally insulates the um- 
bra, whereas heat exchange along the magnetopause is al- 
lowed in subphotospheric layers as long as the slope of the 
magnetopause relative to the vertical exceeds a critical an- 
gle of typically 30°, which corresponds to the bottom of 
the penumbra, z^ p (see Fig. |l|). Since the energy transport 
is more efficient in the quiet sun stratification than in the 
penumbra, the latter is heated. The rate of heating de- 
pends on the transmissivity of the magnetopause, which 
is a free parameter of the model. It is assumed that the 
transmitted heat is distributed instantaneously on hori- 
zontal planes in the penumbra model. Consequently, as 
it is adopted for the umbra and the quiet sun, also the 
penumbral thermal stratification varies only with depth. 
Thus, the surplus brightness of the penumbra in the model 
is due to a finite transmissivity of the magnetopause. The 
tripartite model of JS has been successful in reproduc- 
ing the observed magnetic flux of the spot, the observed 
radial variation of the mean magnetic field strength, the 
observed radii of sunspots, and the Wilson depressions 
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of umbra and penumbra. Moreover, interchange convec- 
tion naturally causes an inhomogeneous penumbra with 
features that are concurrent with many high resolution 
observations (e.g., Schmidt et al. 1992; Title et al. 1993; 
Lites et al. |1993| ). 

The present investigation was started in order to test 
quantitatively the concept of interchange convection. Nu- 
merical investigation of the dynamics of the whole ensem- 
ble of magnetic flux tubes is not possible at present. As 
the natural first step, the discussion presented here is de- 
voted to the evolution of a single magnetic flux tube, which 
is embedded in a tripartite sunspot. Moreover, we assume 
the flux tube to be thin. The equations that govern the dy- 
namics of a thin magnetic flux tube are described in the 
next section. We also describe there briefly the numeri- 
cal realization of the model. Sect. 3 describes the model 
set-up, i.e. the initial conditions and the boundary con- 
ditions, and the resulting evolution of a thin magnetic 
flux tube. In Sect. 4, we compare the results of numeri- 
cal simulations with the observed photospheric features of 
penumbrae, and argue that our model can reproduce the 
penumbral fine structure and the Evershed effect. Also, 
the concept of interchange convection is discussed on the 
basis of these new results. In the last section we summa- 
rize the results emphasizing the observational signatures 
that are suggested by our model. 

2. The model 

We start by writing down the full 3-dimensional system of 
equations for ideal magnetohydrodynamics (MHD) (see, 
e.g., Priest 1982). In ideal MHD, the electrical conductiv- 
ity is assumed to be infinite, leading to frozen-in-magnetic 
fields. Using the substantial (convective) derivative, 



d d 



(1) 



and cgs-units, the equations for continuity, induction, e- 
quation of motion, entropy equation, divergence free mag- 
netic field, and the equation of state read: 



dp 
dt 

B 



d 



-/A7-v , 
V x (v x B), 



dt 

VB = , 

dv 1 
p— = -Vp + gp+ — (V x B) x B, , 



dt 

T dS (AS 



47T 



rati 
4cr 



p 3c 



(2) 

(3) 
(4) 
(5) 

(6) 
(7) 



Here, p, T, p, B, v, and g have their usual mean- 
ing. S denotes entropy per unit mass, 1Z the gas con- 



stant, p = p(p,T) the mean molecular weight. The sec- 
ond term on the right hand side (RHS) in Eq. ^ de- 
scribes the radiation pressure, with a being the Stefan- 
Boltzmann-constant, and c the speed of light. By tak- 
ing advantage of the thin flux tube approximation (De- 
fouw |1976| Spruit |l981a[ |1981bfc Ferriz Mas & Schussler 



1989) this 3-dimensional problem can be reduced to a 1- 



dimensional problem. Doing so, we closely follow the ap- 



proach of Moreno Insertis (1986, hereafter MI). Further 



more, adopting the tripartite model of JS as a background, 
we can describe the dynamics of a thin flux tube evolving 
in a 2-dimensional plane (x,z). 

2.1. Walen equation and equation of motion 

In a thin magnetic flux tube physical variables vary only 
along the tube. Thus, the tube constitutes a curve in a 
plane, and it is appropriate to use a Lagrange represen- 
tation. As the independent variable, we have chosen the 
integrated mass, a, along the tube. 

The position of a mass element is specified by the vec- 
tor: x(a,t) — (x(a,t), z(a,t)). Then the velocity and tan- 
gent vectors are defined as follows: 



v(a,t) = 
l(a,t) = 



dx(a, t) 



dt 

dx(a, t) 



Da 



(8) 
(9) 



Note, that d/dt(...)\ a is identical with the substantial de- 
rivative, d/di, of Eq. (Q). The arc length, ds, of a mass 
element, da, is given by ds = I ■ da, with I = 

Within the thin flux tube approximation, the magnetic 
field vector in the tube is always parallel to the tube, 



B = Bt 



(10) 



with the tangent unit vector, i = 1/1. The magnetic flux, 
0, is conserved along the tube, i.e., 



BA = constant , 



(11) 



with A being the cross section of the tube. With the help 
of Eqs. (||), (|l7j|), and ( |il"| ) one immediately gets, 



— =(j)l. 
P 



(12) 



Finally, one assumes that the total pressure equilibrium of 
the tube with its surroundings is achieved instantaneously: 



P 



£2 
8tt 



= Pb 



8tt 



(13) 



The index b denotes background variables. 

The continuity and the induction equations combined 
together yield the Walen equation which can be written 
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A 



A 
t 



with d denoting the diameter of the tube. 

It is appropriate to decompose the equation of motion 
into a component parallel to the tube, and a component 
perpendicular to the tube. For the parallel component, 
only the first two terms of Eq. (|^) are relevant. Hence, 



A 



^ dv /dv 
' df~ ~~ \~dt 



1 dp , 
pi oa 



(19) 



where 1? denotes the angle between the vertical and the 
tangent vector, t (see Fig. |^). The perpendicular compo- 
nent is obtained by multiplying the final form of Eq. ( |l7| ) 
with the normal unit vector n: 




Fig. 2. Geometric variables. The motion of a curved 
tube within a 2-dimcnsional plane can be characterized 
by the position vector x(a, t). The independent variable a 
denotes the integrated mass along the tube, t denotes the 
tangent unit vector, n the normal unit vector, and $ the 
angle between the tangent vector and the vertical. 

with a use of relations given by Eqs. (^|) and ( |i~2| ) in the 
following form: 



— = (I ■ V)v = — 

dt y ' da 



(14) 



The third term of the equation of motion (||) can be 
simplified using the vector identities: 



(V x B) x B = B 2 nn 



Ud_ 
2Ws 



B 2 



-VB 2 



(15) 



where fi denotes the normal unit vector as depicted in Fig. 
H, and the curvature k is defined as 



nn = 



dt 

ds 



(16) 



Then, using the condition of the instantaneous pres- 
sure equilibrium (Eq. |l3|), the equation of motion can be 
rewritten: 



dv 

dt =g 



B 2 k _ t dB 2 



-n 



47T/9 8np ds 



8np p 



(17) 



The additional term on the RHS, i*b / p, takes into account 
the fluid flow around the tube. i*b denotes the drag force . 
We follow the work of MI (see also Caligari et al. 1995) 
and use 



AF D = -pb-^v±\v±\ 



(18) 



dv 
~dJ 



dv 

dt 



g sin j? ( 1 — 



Pb 



B n 

Anp 8np 



VBl 



Fd 
P 



(20) 



Equations @, ©, and © deter mine the time evo- 
lution of the tangent vector l(a, t) — (l x (a, t), l z (a, t)) and 
the velocity v(a,t) = (v x (a,t),v z (a,t)). Here, v x and v z 
are obtained by projecting the velocity component accord- 
ing to i9 = a,Tctaiil x /l z . These equations are identical to 
those derived by MI, except that our tube is embedded in a 
magnetic background. Therefore, our perpendicular com- 
ponent of the equation of motion contains an additional 
force (3 rd term on RHS of Eq. ©), which results from 
the fact, that the gradient of the magnetic pressure in the 
background exerts a force on the tube. In the magneto- 
hydrostatic background itself this force is canceled by the 
magnetic curvature (2 nd term on RHS of Eq. (p0[)). An- 
other subtle difference lies in Eq. (p"3|): in our case, the 
external pressure is given by the sum of the background 
gas pressure and the magnetic pressure of the background. 

2.2. Entropy equation 

The thermodynamic state of a mass element is character- 
ized by p(a, t), p{a, t), T(a, t), and B(a, t). These variables 
are determined by three algebraic equations and one dif- 
ferential equation, i.e. Eqs. (0), (|i~2|), (|l~3|), and (||). 

The thermodynamic parameters of the flux tube are 
calculated in exactly the same way as the thermodynamics 
in the background (see Jahn 198S for details) including 
the effects of partial ionization on 7 = c p /cv and on the 
mean molecular weight, p = p(p, T) (note that 7 and p are 
not assumed to be constant). The chemical composition 
is assumed according to Stix ( 1989| ): X — 0.7331, Z = 
0.01669. 

Using the first law of thermodynamics we can write 
the entropy change in the following form: 



J_dS 

Cy dt 



dhxT 
dt 



(7-1) 



Xp _ dlnp 
XT dt 



(21) 
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with xt '■- 
and x P - 



d lnp 



dlnT 
dhip 



9 hip 



(22) 
(23) 



This expression is substituted into the entropy equation 
(|J). After eliminating the total derivatives dlnpldt and 
din T/dt, that result from differentiating Eqs. (Q), (|l2|), 
and (U3), one finally obtains 



-75 ~ 7Xp 



2 

c v \dt / 



din, 
~df 



rad 



ldP b 2 dln^ 
p~dT + P dt 



(24) 



Note that the background is hydrostatic, dP^/dt = 0, so 
that dPb/di = (v • V)Pb, i-e., the temporal change of P b 
is only due to the motion of a Lagrange point in the 2D 
plane. 

2.2.1. Heat exchange by radiation 

Since radiation is the dominant physical process in the so- 
lar photosphere, radiative heat exchange of the tube with 
the surroundings has to be taken into account. We do 
this by adopting the relaxation time approach by Spiegel 
( |1957| , al so use d by Degenhardt |l99l| and Thomas & Mon- 
tesinos ( 1993 )). For a homogeneous atmosphere in local 
thermodynamic equilibrium, the inverse damping factor 
of a temperature perturbation can be identified with a 
radiative relaxation time. 



16kaT 3 



arccot(r)] 



(-1) 



(25) 



where k is the Rosseland mean opacity (Huebner et al. 



1977). For small differences between the two temperatures, 



the entropy change by radiation is then given by, 



dS 
~dt 



= -i2-(T b -T) 

tad r rad-t 



(26) 



2.3. Numerical treatment 



The evolution of a thin magnetic flux tube is given by 5 
differential equations (two components of Eq. ( fL4|) , (|l9|), 
(2C), and (p^)) and 3 algebraic equations ((0), (|l2j), and 
(13)), which are solved for the 8 unknown variables Z(a, t), 
v(a, t), p{a, t), p(a, t), T(a, t), and B(a, t). The tube's cur- 
vature, k, is given by 



1 ( dl x _ di z 

P'\ z da x da 



(27) 



It is appropriate to use a staggered grid (see MI), in 
which x and v are discretized on the boundaries of a cell, 
and the thermodynamical state, i.e. p, T, p, B, and 1 in the 



Table 1. Parameters of background model (cf. JS) 



Parameter 



Symbol Value 



total magnetic flux 


$tot 


1 • 10 22 Mx 




$ 


n 75 <j>t * 

VJ. 1 >j ^tot 


Wilson depression of umbra * 


W u 


-470 km 


Wilson depression of penumbra 




-150 km 


radius of penumbra at z = W pu 


-Rpu 


13 500 km 


radius of umbra at z = W u 


-fill 


5000 km 


transmissivity 


C 


0.7 


lower boundary of model 


Zbot 


-20 000 km 


upper boundary of model 


Ztop 


800 km 


bottom of penumbra 


Zbp 


-6 000 km 


B at the bottom of model 


Bbot 


15 kG 


B at the center in photosphere 


B ccn 


2700 G 



' z = km corresponds to the photosphere (r = 2/3) in 
the quiet sun and decreases inwards. 



center of the cells. The time integration is done explicitly 
in two steps. First, the quantity / is calculated for an in- 
termediate step, f n+1 ' 2 = f n + {%) n -dt/2. Then, the full 
time step is calculated, f n+1 = f n + (^ ) n+1/2 ■ dt. The 
time step used adapts itself, securing that the Courant- 
Friedrichs-Lewy-criterion for stability is satisfied for each 
cell, dt < min(ds/t) s , ds/vA). Here, ds denotes the ar- 
clength of one cell, v s = \Jxp c pl c v • \pp~l P the sound 
velocity, and «a = B /s/iTrp the Alfven velocity. Simul- 
taneously, the time step is forced to be smaller than the 
radiation relaxation time, t Ta ,d- Once, l n+1 , v n+1 , and p n+1 
are known, B n+1 is determined by Eq. (|l^) and the posi- 
tion vector, x n+1 , that fixes the background variables, is 
obtained by integrating the velocity in time. With B n+1 , 
p b i+1 , and P™ +1 , Eq. @ yields the gas pressure, p n+1 . 
The Saha equations determine the ionization fractions of 
hydrogen and helium, and thereby the mean molecular 
weight fi and the specific heat ratio 7. The proper pair 
of ix n+1 and T n+1 is found by iterating the equation of 
state, given by Eq. (R). The Saha equations and the itera- 
tion scheme are described in more detail in Schlichenmaier 



( |1997| ) 



2.4- Embedding the tube in a tripartite model 

To follow the evolution of the flux tube in the penumbra 
we embed it in a tripartite sunspot model. Here, the idea 
is that we take a bundle of magnetic field lines from the 
sunspot model. These field lines are initially in magneto- 
hydrostatic equilibrium and form a physical entity, being 
described by a thin magnetic flux tube in ideal MHD. The 
parameters of the particular tripartite model we used are 
listed in Table 0, for an illustration see Fig. 1 in JS. 

The thermodynamic stratifications for the umbra, the 
penumbra, and the quiet sun are prescribed by the tri- 
partite model. Within these stratifications the pressure, 
Ph{z), the density, pb(z), and the temperature, T^(z) de- 
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Temperature T [K] 



Gas pressure p [dyu cm 




geometrical depth z in Mm 



Fig. 3. Thermodynamic variables of the back- 
ground. A sunspot is approximated by three different 
stratifications, i.e., the umbra (thick line), the penumbra 
(regular line), and the quiet sun (thin line). 




«dial distance x in Mm 



Fig. 4. The background magnetic pressure. B^/8n, 
is shown as a logarithmic surface plot versus the radial dis- 
tance, x, and depth, z. For the quiet sun, we set log(O) = 1. 



pend only on depth. Level z = Okm corresponds to 
r = 2/3 level of the quiet sun and decreases downwards. 
Fig. H shows the dependences of Tb, pb and pb on z for 
the different stratifications in the vicinity of the photo- 
sphere. The thermodynamic variables inside the tube, at 
the center of each cell, are obtained with a use of logarith- 
mic splines defined on the background discretization. The 
tube lying along the magnetopause is illuminated by the 
quiet sun on one side. This effect is taken into account in 
the radiative heat exchange. 



The background magnetic field strength, Bb(x, z), de- 
pends on the radial distance, x, and on depth, z. The 
radial distance x is defined to be zero at the center of the 
sunspot model. Bb(x,z) for a certain position is deter- 
mined using a bilinear interpolation. The magnetic pres- 
sure, B^/8tt, of the used background model is shown in 
Fig. | as a surface plot. One can clearly see the jump 
across the magnetopause from finite field strength in the 
penumbra to a vanishing field in the quiet sun. The current 
sheets, i.e. the magnetopause and the peripatopause, con- 
stitute discontinuities for the magnetic pressure, B^/8n, 
the gas pressure, pb, the temperature, Tb, and the den- 
sity, pb- But note, that when crossing the current sheets 
horizontally, the total pressure (gas plus magnetic) is con- 
stant. 

3. The evolution of a thin magnetic flux tube 
embedded in a sunspot 

An ensemble of magnetic flux tubes is assumed to partic- 
ipate in interchange convection. As a first step we restrict 
the study to simulations of the evolution of a single flux 
tube. In this contribution, we present the results obtained 
for one specific flux tube and concentrate mainly on its 
photospheric behavior. 

3.1. Initial configuration 

The results presented in this section have been obtained 
for the tube with a magnetic flux <fr = 2 • 10 16 Mx. That 
corresponds to a diameter of d = 2 • \J (f)/Bir = 50 km, for 
a magnetic field strength equal to 1000 G. Initially, the 
tube lies along the magnetopause. Hydrostatic equilibrium 
along the tube is satisfied since the distribution of the 
density and temperature is identical with the penumbral 
stratification: 

t = : p(x, z) = PKpcnumbra^) 
T{x, z) — T'b ? penumbra('^) 
p{x^ z) — £*b, penumbra^) 

Magnetostatic equilibrium, i.e., the cancellation of the sec- 
ond and third term on the RHS of Eq. (20), is an intrin- 
sic property of the tripartite model, and can be realized 
numerically with an accuracy which is sufficient, but lim- 
ited by the interpolation procedures. The magnetic field 
strength is given by the total pressure equilibrium across 
the magnetopause, 



t = 0: B(x,z) = V87T-y/pb, quiet sun(z) ~ Pb,pcnumbra(^) • 

The model-tube extends down to a depth of z = —15 
Mm at a radial distance of x = 4.8 Mm. The upper end 
is placed at a radial distance of x — 24 Mm and at a 
height of z = 0.3 Mm above the photosphere (r = 2/3) 
of the quiet sun. The intersection of the tube with the 
photosphere coincides initially with the outer edge of the 
penumbra model. 
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3.1.1. Boundary conditions 

The lower end of the tube, which is placed at the local 
position of the magnetopause at z = —15 Mm deep in 
the convection zone, is held fixed, i.e. the mass element at 
the lower boundary has zero velocity for all times. At the 
upper boundary, we apply a free boundary condition: The 
time derivative of the tangent vector, dl/di, is extrapo- 
lated and the total force at the last grid point equals the 
total force at the last but one grid point. In case of an 
outflow, mass elements reaching x — 27 Mm are cut off. 
That implies that the number of grid elements gradually 
decreases, if an outflow is present. 

3.1.2. Numerical tests 

The results presented in this section were obtained using 
1033 grid points, being distributed equidistantly in mass 



for z > 



itch 



T Mm and equidistantly in arclength 



below that depth. It has been tested that the evolution is 
the same for finer grids and does not depend on t he par - 
ticular choice of the depth, z S witch (Schlichenmaier 1997 ). 

When the radiative heat exchange is switched off, en- 
tropy becomes a constant of motion. Since the entropy of 
the system is given by its thermodynamic properties, the 
conservation of entropy per discretized cell can be checked 
independently from the code (Schlichenmaier 1997 ). For 
this run, the code preserves entropy to a relative accuracy 
of less than 10~ 4 . 

3.2. Time series 

Fig. H shows the time series of our simulation. Each im- 
age shows the shape of the tube. The diameter of the 
tube has been magnified by a factor of 12, for better visi- 
bility. The gray-coding represents the velocity field along 
the tube. The background-shading represents zero veloc- 
ity, while darker coding corresponds to positive velocities 
which point upwards and outwards. Images are plotted 
every 600 s. 

3.2.1. The onset of evolution caused by radiation 

Even though initially the tube is in magnetohydrostatic 
equilibrium, it will start to rise because of the radiative 
exchange of heat. Part of the tube is illuminated by the 
underlying hotter quiet sun so that it heats up and ex- 
pands relatively fast. That leads to a decrease of density 
inside the tube as compared to the outside, and the buoy- 
ancy (first term on RHS of Eq. (p0[)) causes the tube to 
rise. The region where the ascending motion appears first 
is determined by an interplay between the radiative relax- 
ation time, i ra d, the temperature difference, T{z) — T\,(z), 
and the superadiabaticity, 5 — V — V a d. Since radiation is 
most effective near the photosphere, the radiative relax- 
ation time (see Eq. (|2E|)) is the shortest there. Superadia- 
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Fig. 5. Time series of the evolution of a flux tube that 
contains <f> — 2 • 10 16 Mx, and initially lies along the mag- 
netopause. Starting just below the photosphere the tube 
rises, thereby causing the footpoint of the tube to migrate 
radially towards the umbra. While it rises, a flow along 
the tube develops, being indicated by the gray-coding of 
the tube. The background-gray corresponds to km s . 
The flow is directed upwards and outwards. At the outer 
edge of the penumbra the flow velocity becomes as high 
as 14 km s _1 . 

baticity, being the cause of convective instability, is largest 
in the surface layers - within the first few hundred kilo- 
meters below the photosphere. The temperature difference 
between the quiet sun and the umbra, being a measure for 
the rate of heating, is largest near the photospheric level of 
the penumbra, i.e., at z — —150 km. Therefore, the tube 
gains buoyancy at the fastest rate just below the photo- 
sphere and there the tube will start to rise first. That can 
be seen in the second image of Fig. HJa. 

3.2.2. Initial stages of the evolution 

Below the photosphere, the stratification is convectively 
unstable (superadiabatic) and the tube's rise accelerates. 
Above the photosphere, the stratification is convectively 
stable, and moreover, the tube's density increases due 
to radiative losses. Therefore, the tube's rise decelerates 
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rapidly in the atmosphere. Approximately 100 km above 
the photosphere, a new equilibrium is found for the tube, 
in which the downward-acting buoyancy force is neutral- 
ized by the upward-acting gradient of the background 
magnetic pressure (third term of RHS in Eq. fl2C|)). Simul- 
taneously, the footpoint of the tube (i.e. the intersection of 
the tube with the photosphere at r = 2/3), moves radially 
inwards toward the umbra. 



3.2.3. Build-up of surplus gas pressure inside the rising 
tube 

A positive flow velocity along the tube develops as the 
tube rises. It reaches a maximum of 14 km s _1 at the 
outer edge of the penumbra at later stages of the evo- 
lution. The onset of this longitudinal flow results from 
an enhanced gas pressure gradient along the tube, which 
builds up locally at the place where the tube rises into 
the photosphere. In order to understand this phenomenon, 
one has to realize that the tube finds itself (locally) in a 
strongly superadiabatic environment characterized by a 
very small gas pressure scale height (of the order of 100 
km) which is much smaller than the scale height of the 
magnetic pressure (approximately 3 000 km). Note also, 
that below the photosphere the gas pressure dominates 
the magnetic pressure, i.e. (3 > 1. 

The tube stays in total (gas + magnetic) pressure equi- 
librium as it rises (cf. Eq. (13)) and it must expand, be- 
cause of the decreasing pressure in the background. Since 
the decrease of the total background pressure is dominated 
by the gas pressure, the magnetic pressure, B 2 /8ir, inside 
the tube gets smaller than the background magnetic pres- 
sure, B 2 /8n. In consequence, the gas pressure inside, p, 
must build up relative to the background gas pressure, 
Pb- Note, that the background gas pressure is constant 
horizontally within the penumbra. 



Since the part of the tube lying above the photo- 
sphere does not rise, it becomes almost horizontal be- 
tween the footpoint, at Xf p , and the outer edge of the 
penumbra, at x oc , i.e. ph(x[ p ) ~ Pb(x oc ). Near the outer 
edge of the penumbra the gas pressure inside the tube 
does not change, p(x oc ) = Ph(xoe), whereas the gas pres- 
sure is greatly enhanced in the vicinity of the footpoint, 
p(x[ p ) > pt^Xfp), since there the tube has risen through 
the strongly superadiabatic region. The resulting gas pres- 
sure gradient accelerates a horizontal outflow along the 
tube which is accompanied by an upflow in subphoto- 
spheric layers. This upflow is maintained by the superadia- 
batic background stratification. The gas pressure gradient 
in the horizontal part of the tube is sustained by radia- 
tive losses which decrease the internal energy (and the gas 
pressure) of the matter flowing along the tube. 



3.2.4. Centrifugal force and magnetic tension 

The migration of the footpoint is governed by three forces 
acting perpendicular to the tube (the gradient of the back- 
ground magnetic field can be neglected here). First, at the 
turning point, i.e. slightly above the footpoint, where the 
tube bends horizontally, magnetic curvature force is pos- 
itive (cf. Fig. ||) and it decelerates the migration of the 
footpoint. Second, the longitudinal flow exerts a centrifu- 
gal force at the turning point. The centrifugal force term 
stems from the advection term in Eq. ([j]), v ■ V: 



Vv 



2 ds 



v 2 Kn 



The first term gives the inertia force along the tube, and 
the second term is the centrifugal force, which is the inertia 
force perpendicular to the tube. 

Third, the buoyancy being a vertical force that acts 
perpendicular to the tube dominates the rise as long as 
the tube is highly inclined with respect to the vertical. As 
the time series shows, the inward migration of the foot- 
point ceases at t = 3 600 s. Here, the buoyancy becomes 
negligible, because of the almost vertical orientation of the 
subphotospheric part of the tube. Thus, in later stages of 
the evolution, the perpendicular component of the equa- 
tion of motion near the turning point is governed only by 
the centrifugal and the magnetic curvature forces. It is in- 
teresting to note, that the centrifugal force neutralizes the 
magnetic curvature force when the flow velocity equals the 
Alfven velocity. 

3.2.5. Formation of a shock front above the photosphere 

The time series shows a shock front that develops beyond 
the outer edge of the penumbra. It forms at x = 15 Mm 
after 1 800 s and migrates outwards with a velocity of 
roughly 3 km s . While it migrates outwards and up- 
wards, the shock front accelerates and reaches ss 6 km s _1 
at the outer (upper) boundary. Based on simple physical 
arguments, the mere existence of this shock wave can be 
explained as follows. Assume you have a stationary flow 
upwards along the tube. Due to continuity, the flow veloc- 
ity, vu , would then increase according to uii ~ 1/ p, whereas 
the Alfven velocity va scales as va ~ B/y/p. Now, p de- 
creases exponentially and, not so important, B decreases 
linearly with height. Thus, vi\ being smaller than va at 
some initial height, increases more rapidly than va- At 
some point, vu has to become superalfvenic. In a magnetic 
flux tube the critical flow velocity along the tube is given 
by the tube speed, v t — va- v s /^/va + v s (e.g., Ferriz Mas 
1988), implying that a superalfvenic flow is supercritical 
at the same time. Having a supercritical flow that collides 
with plasma which is at rest initially has to result in a 
shock front. As this shock-front accelerates the plasma, it 
migrates outwards and leaves the computational domain 
after 4 200 s. 
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In the present contribution we do not make any at- 
tempt to realistically describe shock waves. But it is clear 
that the shock does not influence the dynamics happen- 
ing near the pcnumbral photosphere for several reasons. 
First, the shock front is a consequence of the acceleration 
of plasma within the penumbra, rather than the cause of 
the longitudinal flow. Second, the density decreases expo- 
nentially with height, and therefore the inertia within the 
photospheric penumbra is much to large to be affected by 
the dynamics happening hundreds of kilometers above the 
photosphere. And after all, the flow velocity upstream of 
the shock is supercritical, i.e., the penumbral photosphere 
doesn't know about the shock front. 

The flow is decelerated within the shock and the gas 
pressure increases there. That leads to a local expansion of 
the tube at the shock front. Such an expansion should be 
counteracted by magnetic tension forces inside the tube. 
However, within the thin flux tube approximation such 
forces are neglected and nothing prevents the tube from 
infinite expansion at the shock front. Therefore, we have 
included the magnetic tension forces inside the tube ficti- 
tiously by changing the background pressure, Pt,: At the 
location of expansion we enhance the background pres- 
sure. The stronger the expansion, the larger the back- 
ground pressure. 



3.2.6. Final stages of the simulated evolution 



As discussed in Sect. 3.2.3, plasma within the tube ex- 
pands as it rises through the convection zone. Thereby, the 
magnetic field strength, B, decreases. The magnetic field 
strength of plasma reaching the footpoint is the smaller, 
the larger is the height difference it traveled through the 
convection zone. That implies, that va = B / ^/AWp at the 
footpoint decreases continuously, while the flow velocity 
stays constant in time at the footpoint. Hence, at some 
point (in our simulation at t s» 7000 s) the flow becomes 
superalfvenic and the centrifugal force exceeds the mag- 
netic tension force at the turning point, i.e. the flow over- 
shoots the turning point. Due to radiative losses in the 
atmosphere, the plasma that overshoots gets denser and 
decelerates. 

The amplitude of the overshoot gradually increases in 
time as Vh/va increases. Once the flux tube dives back 
into the convectively unstable subphotospheric penumbra, 
anti-buoyancy drags the tube down, and magnetic tension 
is too weak to prevent the tube from sinking. Fig. |^b shows 
a snapshot of the tube in the vicinity of the turning point 
after 9 300 s. The inertia of the upflow at the footpoint 
prevents the turning point from sinking. The part of the 
tube that dives back beneath the photosphere sinks down. 
Since the radii of curvature at the turning point (x « 8.4 
Mm) and at x w 10.2 Mm become very small, the thin flux 
tube approximation is no longer valid. For that reason, we 
stop the simulation here. 
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Fig. 6. Two snapshots of the flux tube and the 
temperature profile. As the Alfven velocity increases 
at the footpoint due to the continuous upflow along the 
tube, the centrifugal force at the turning point (x w 8.4 
Mm) exceeds the magnetic tension, leading to overshoot- 
ing, a) shows the phase in which the overshoot is convec- 
tively stable at t = 7 500 s, while b) visualizes the instable 
phase at t = 9 300 s, where the tube dives back into the 
subphotospheric stratification. The gray coding refers to 
temperature according to the shading bar. 



4. Observational consequences 

In order to discuss the observational consequences of our 
model, we concentrate on the evolution of the tube in and 
above the photosphere, since only here is the tube ob- 
servable. Fig. Q shows the proper inward velocity of the 
migrating footpoint as a function of the radial distance 
x from sunspot center. It illustrates the migration of the 
tube's footpoint which will be compared with the observed 
inward migration of bright penumbral grains. Further, we 
use Fig. |^ to discuss the features of our simulations in 
the context of the observed penumbral fine structure such 
as bright and dark filaments and the Evershed effect. It 
shows an intermediate stage of evolution of the tube near 
the surface layers. From left to right, one can see the um- 
bra, the penumbra, and the quiet sun being separated by 
the peripatopause and the magnetopause, respectively. In 
Fig. ||a, ||c, and ||d the gray coding represents the temper- 
ature variation, the gas pressure variation, and the varia- 
tion of the magnetic field strength, respectively. Note, that 
in each background stratification the temperature and the 
pressure vary only with depth, whereas the magnetic field 
is two-dimensional. Fig. ||b shows the optical thickness of 
the tube. Here, the subphotospheric background is illus- 
trated as an optically thick regime and the atmosphere as 
an optically thin regime. For better visibility the diame- 
ter of the tube in Fig. || has been magnified by a factor 
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that the tube's diameter decreases between the footpoint 
and the point of thermal equilibrium. Thus, our model is 
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Fig. 7. Proper motion of the footpoint. The appar- 
ent inward velocity of the migrating footpoint is plotted 
versus the radial distance x from sunspot center. The plus 
signs are equally spaced in time (At — 300s), tracing the 
evolution between t = 300 s (x = 13 300 km) and t — 7 200 
s (x = 8 300 km). 



6. The arrows inside the tube indicate longitudinal flow 
velocities. Just below the footpoint a velocity of 3 km s _1 
is present, beyond the footpoint the velocity increases to 
6 km s _1 , and at the outer part of the penumbra near 
x = 12 Mm the flow velocity reaches 13 km s _1 . 



4-1- Penumbral grains and bright filaments 

The footpoint of the flux tube can be identified with a 
bright penumbral grain, because the upflow along the tube 
brings hot plasma to the surface which appears bright 
against the darker background. As demonstrated in the 
previous section, the tube's footpoint migrates from the 
outer edge of the penumbra inward towards the umbra. In 
Fig. the data points (plus signs) are equally spaced in 
time (At = 300 s), illustrating that the footpoint starts 
off at x = 13 300 km with a proper velocity exceeding 2 
km s" 1 . After 1 200 s (4 th data point from the right) the 
migration of the footpoint decelerates from 1.5 km s _1 to 
« 0.1 km s" 1 at t = 7 200 s. These values arc consistent 
with the proper motion of penumbral grains observed by 
Muller ( |l973b| ). 

At the footpoint, the tube's temperature and gas pres- 
sure are higher than the corresponding background val- 
ues. As the plasma flows along, the tube loses internal 
energy by radiation. Thereby it cools and the gas pres- 
sure diminishes. At x ~ 10 Mm it reaches temperature 
equilibrium with the surroundings. As the gas pressure 
diminishes, the magnetic field strength increases accord- 
ing to Eq. (|T^), and flux conservation, Eq. (|ll|), implies 



concurrent with observations from Muller (1973a, b) and 
Tritschler et al. (1997) which show that penumbral grains 
consist of a bright coma and a somewhat dimmer thinner 
tail which is radially elongated and points away from the 
umbra. Our model suggests that bright filaments are the 
long tails of penumbral grains. However, for lower spatial 
resolution, a few radially aligned penumbral grains might 
also be visible and interpreted as one bright filament. We 
want to mention, that the length of the bright tail de- 
pends on the amount of magnetic flux of the tube: More 
magnetic flux implies a larger diameter and thus a larger 
optical depth, a longer radiative relaxation time (see Eq. 
25])), and hence a longer bright tail. 



4-2. Dark filaments 

In view of Figs. ||a and ||b, one can see that the tube is 
optically thick as long as it is hotter than the surround- 
ings. This is due to the temperature dependence of the 
H~ opacity which is the dominant opacity near the pho- 
tosphere: k oc T 9 . Another important feature was already 
mentioned in Sect. 3.2.2: The tube is elevated above the 
photosphere. That means that bright filaments are opti- 
cally thick structures overlying a darker background and 
surrounded by an optically thin atmosphere. Therefore, we 
propose that dark filaments might not exist per se. They 
rather appear as a consequence of spacing between r adiall y 
elongated bright filaments (Schlichenmaier et al. 1998 ) 



This i s consistent with the statement of Muller ( 1973a . 



1973b) that within the penumbra bright features show up 



against a dark background. 
4-- 3. Evershed effect 

Our model also offers a consistent explanation for the Ev- 
ershed effect: Between the footpoint, at x « 10.5 Mm, 
and the outer edge of the penumbra, the matter inside 
the tube becomes gradually more transparent. The opti- 
cal thickness of the tube decreases to t m 10 _1 (see Fig. 
||b) and a line of sight that crosses the tube reaches op- 
tical depth r = 2/3 in the photosphere of the model at 
z — —150 km, i.e. below the tube. In other words, after 
the plasma has cooled reaching temperature equilibrium 
with the surroundings, the tube becomes transparent. The 
tube constitutes a flow channel that is elevated approxi- 
mately 100 km above the photosphere and confines an 
outward flow reaching velocities up to 13 km s . In the 
following we discuss the influence of such a flow channel 
on the line profiles of photospheric absorption lines. 

Generally speaking, a flow channel Doppler-shifts a 
line profile that is formed inside the flow channel. In order 
to understand the line asymmetry which is a characteris- 
tic feature of the Evershed effect, one should recall that 
the contribution function of the line core has its maxi- 
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Fig. 8. a-d. Snapshot of the evolution at t = 5 400 s. From left to right one sees the umbra, the penumbra, and 
the quiet sun. The tube that has risen from the magnetopause lies horizontally above the photosphere, being elevated 
some 100 km. The tube's diameter is magnified by a factor 6 for better visibility. As indicated by the legends, the 
gray coding represents physical variables along the tube and in the background: a Temperature, b Optical thickness 
of the tube, c Gas pressure, d Magnetic field strength. The arrows in the tube visualize longitudinal velocities in the 
tube. In the lower left corner of each plot an arrow with a length that corresponds to 5 km s _1 is shown. 



Schlichenmaier, Jahn, Schmidt: Magnetic flux tubes in sunspots 



13 



contribution functions j ^ 



line core 



z in km 



line of sight 




photosphere 




line wing , 



Fig. 9. Explanation of line asymmetries. A thin flow 
channel only affects parts of a typical photospheric ab- 
sorption line. As it is drawn in this figure, the flow chan- 
nel Doppler-shifts a fraction of the contribution function 
for the line wing, but does not influence the contribution 
function of the line core, leaving the line core unshifted. 
Hence, a thin flow channel can naturally explain the ob- 
served line asymmetry which is a characteristic feature of 
the Evershed effect. 



mum higher up in the atmosphere than the contribution 
function of the line wing, as it is depicted for a typical pho- 
tospheric absorption line in Fig. ^. Since the flow channel 
might be substantially thinner than the width of the cor- 
responding contribution function, only distinct parts of a 
line profile might be shifted. Therefore, a line asymmetry 
would be caused naturally by a thin flow channel. A line 
asymmetry which is such that the line wing is more shifted 
than the line core (see Fig. lb of Degenhardt 1993 ) can be 
explained by a flow channel that is sufficiently thin and is 
elevated onl y slig htly above the photosphere (Schlichen- 
maier et al. 1998): Such a flow channel affects the contri- 



bution function of the line wing, but does not influence 
the contribution function of the line core which originates 
in still higher layers. 

The diameter of the tube in our simulations is less than 
50 km, depending on the strength of the magnetic field. 
This has to be taken into account, when shifted asym- 
metric line profiles are interpreted in order to deduce a 
flow velocity. Wiehr (1995) deduces a flow velocity of > 5 
km s _1 by assuming that the line asymmetry is caused by 
a spatially unresolved line satellite being caused by a flow 
channel. If one further takes into account, that a typical 
contribution function has a full width at half maximum of 
more than 100 km and that our simulated tube has a di- 
ameter of less than 50 km, a local flow velocity exceeding 
10 km s _1 may still be consistent with observation. 



4-4- Outward flow in bright filaments 

In the preceding section we have demonstrated that the 
outward flow in the transparent cooler part of the tube can 
reproduce a line profile that is shifted and asymmetric. 
However, our model also predicts an outward flow in opti- 
cally thick bright filaments. This result is somewhat sur- 
prising, since it is observed that the Evershed flow is corre 
lated with dark filaments (e.g., Beckers & Schroter 19i 
Title et al. |1993|; Wiehr & Degenhardt [1994D and that 



bright filaments are anti-correlated with the Evershed flow 
(Shine et al. 1994). We explain this apparent contradiction 
by the fact that all these measurements were done using 
Fe I lines that have a low excitation potential of less than 4 
eV. Since these lines originate almost exclusively at atmo- 
spheric temperatures of, say, less than 6000 K, we expect 
that a hot, optically thick flow channel has a negligible 
contribution to the intensity profile of these lines. If any- 
thing, a hot flow channel would at most shift the line wing, 
where the line intensity is close to the continuum inten- 
sity. Hence, a hot flow channel may not be detectable with 
absorption lines of low excitation potential. On the other 
hand, absorption lines which arc sensitive to high temper- 
atures, i.e. lines which have a higher excitation potential, 
are expected to originate mainly within a hot flow chan- 
nel. Hence, observing bright filaments, i.e. hot plasma, it 
is more appropriate to use a "hot" absorption line with a 
high excitation potential, such as the CI 538.03 nm line. 
Since this line has an excitation potential of 7.68 eV, its 
mere existence within the penumbra indicates that hot 
plasma is present. In fact, high resolution measurements 
in the CI 538.03 nm line demonstrate the existence of this 
line within the penumbra and furthermore, they show line 
core shifts (Rimmele |1995| , Stanchfield et al. |1997[ ). Fig. 5 
in Stanchfield et al. shows line core shifts corresponding 
to Doppler velocities up to 2kms -1 within the penum- 
bra, being blue-shifted on the center side and red-shifted 
on the limb side. We take this observation as an indication 
for the existence of outward flows in bright filaments. 



Finally, we w ant t o mention that Wiehr ( 1995 ) and 
Balthasar et al. ( 1997 ) find indications for the existence 
of a non-zero velocity in bright filaments, if line profiles 
are interpreted in terms of a largely shifted line satellite 
that stems from dark filaments and a weakly shifted main 
component that stems from bright filaments. 

In order to decide whether a line shift would be pro- 
duced by the hot part of our model tube, one has to 
compute synthetic line profiles. Therefore the tempera- 
ture gradient across the tube is needed, which forms when 
an optically thick, hot tube radiates in a cooler, optically 
thin atmosphere. Thus, for computing synthetic line pro- 
files, it is essential to treat the radiative transport for our 
model tube. 

In summary, our model suggests that there are three 
different structures within a penumbra: a) a bright com- 
ponent, in which the optically thick and hot flow may only 
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be detectable using "hot" absorption lines, b) a dark com- 
ponent, in which a flow channel of small optical thickness 
lies above the dark background, and c) a dark component 
consisting of the dark background with no flux tube above 
it. In the last component one would see a dark filament 
with no Evershed flow present, whereas in case b) one 
would look at a dark filament which contains the flow. 
Hence, the Evershed signal is correlated with some but 
not all dark filaments. 

5. Discussion 

We have presented a numerical simulation of the time de- 
pendent dynamic evolution of a thin magnetic flux tube in- 
side a sunspot penumbra. The results yield an extension of 
the studies reported by Jahn et al. ( 1996| ), and Schlichen- 
maier et al. (1997). Here, we aim at a more detailed under- 



standing of the physical structure of the penumbra. Since 
our model seems to reproduce the observed small-scale 
structure of the penumbra, we take this as an indication 
that it supports the concept of interchange convection, on 
which our results are based. The model suggests: 

1. Penumbral grains are the footpoints of inwards mi- 
grating flux tubes, in which the enhanced brightness 
is caused by a systematic upflow. At the footpoint an 
upflow of 3 km s _1 is present. 

2. A bright filament, when interpreted as the dimmer and 
thinner tail of a penumbral grain, is the consequence of 
a hot plasma outflow that cools down during horizontal 
outward motion. Bright filaments are hotter and opti- 
cally thicker than the surroundings, and are elevated 
relative to the dark background. In bright filaments a 
plasma flow is present having a velocity of typically 6 
km s _1 . 

3. The dark background is partly eclipsed by bright fila- 
ments. Hence, the appearance of dark filaments would 
be the result of spacing between bright filaments. 

4. The Evershed flow occurs with velocities up to 13 
km s _1 in optically thin (r « 10 _1 ) and almost hori- 
zontal channels which are elevated some 100 km above 
the photosphere. 



which continues up into the canopy, as it was observed by 
Solanki et al. fll994| ). 

So far, we have only considered the rise of a tube that 
lies along the magnetopause initially. Clearly, for an un- 
derstanding of interchange convection, we need not only 
study different initial conditions, but also the dynamics 
of sinking magnetic flux tubes. If no sinking tubes were 
present, magnetic flux would be transported continuously 
from the outer to the inner part of the penumbra. In conse- 
quence, the magnetic pressure would increase in the inner 
part of the penumbra. We surmise that at some point, 
the gradient of the magnetic pressure, VB^/87r, becomes 
large enough to suppress the further rise and to initiate 
the sinking of the flux tube. 

We want to stress that the acceleration process that 
causes the outward flow is relate d to the siphon flow mech- 
anism (Meyer & Schmidt 1968), in the sense that a gas 
pressure gradient accelerates the plasma. However, in our 
model the pressure gradient is created locally within the 
penumbra: The flow is initiated by the rise of a magnetic 
flux tube and is driven by the superadiabatic background 
stratification. In the horizontal part of the tube, radiative 
cooling sustains the gas pressure gradient that accelerates 
the plasma from 3 km s" 1 near the footpoint up to 14 
km s _1 near the outer edge of the penumbra. In result, 
our model does not require strong magnetic fields at the 
outer footpoint as in the classical siphon flow scenario. We 
find an acceleration mechanism that is caused naturally in 
a certain phase of the evolution of an emerging magnetic 
flux tube. 

Since our simulations yield upflows along subphoto- 
spheric flux tubes, it is tempting to speculate that the sur- 
plus brightness of the penumbra as compared to the umbra 
is caused by these upflows which constitute a convective 
heat transport. Thus, one would suppose that penumbral 
grains and bright filaments which are fed by these upflows 
cause small fractions of the penumbra to be brighter than 
the umbra. 
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In our model, the magnetic flux tube and the corre- 
sponding plasma flow extends beyond the outer penum- 
bral boundary. At the boundary, the tube follows the mag- 
netopause, which slowly gains height and reaches z — 700 
km at x = 27 Mm. There, the boundary condition al- 
lows for a steady outflow of mass. Thus, our model does 
not describe magnetic flux tubes that dive back beneath 
the photosphere near the outer edge of the penumbra, as 



recently reported by Westendorp Plaza et al. (1997) and 
Stanchfield et al. (1997). Modeling dynamical flux tubes 
that bend back down near the outer edge of the penumbra 
requires different upper and lower boundary conditions. 
Instead, our model describes the part of the Evershed flow 
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